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Abstract: Porosity is a major defect in cast aluminum alloys affecting in particular, the fatigue strength. In this work, 
fractal analysis of the pores in the microstructure of cast aluminum alloys were done to provide information linking their 
composition and processing to the crack initiation. 

Two cast alloys: Al - 20%wtSi and Al - 20%wtCu, were used and the types of pores were studied in the as-cast 
samples. The shapes of the pores along with the percentage porosity in each of the micro structures were also evaluated. The 
Multi-Stage Random Sampling (MRS) and Spatial Point Pattern (SPP) Methods were used to determine the distribution of 
the pores and the point of crack initiation. 

Fractal analysis showed that all the pores were shrinkage pores with fi<0.3 and D approaching 2. The MRS and 
SPP methods revealed that crack initiation for eventual failure will start in the worst pore found in the lower right region of 
as-cast Al-20%wtSi because it has the value ofD— 1.2846 and lowest sphericity =0.0011. 

This work has shown the effectiveness of using the fractal analysis, the MRS and the SPP methods for the 
characterization of the pores in the microstructure of cast aluminium alloys. 
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I. Introduction 

For alloys and composite materials containing regular microstructures, a prediction of mechanical properties can 
be made by a quantitative measurement of features such as grain size, particle size, and spacing. This however is not the case 
where an irregular microstructure is involved because of the difficulty in the numerical characterization of the structure. For 
irregular microstructures, the application of fractal geometry offers a method by which both the individual particle shape and 
the mode of the distribution of the particles can be fully described in a numerical manner (1). 

A Study carried out to evaluate the dependence of surface fractal dimension of Al 2 0 3 -Si0 2 composite membranes 
on chemical composition and sintering temperature (2) has shown that all composite membranes have rough surfaces with 
fractal dimensions ranging from 2 to 3. As Si0 2 content increases and sintering temperature rises from 200 to 600°C, the 
fractal dimension increases. However, at 800°C, the surface fractal dimension decreases. The computer aided stereo 
matching method on metals and ceramics was used to reconstruct the three-dimensional images of fracture surfaces formed 
by different mechanisms (3). An extensive experimental investigation of the scaling properties of fracture surfaces in 
heterogeneous materials was done (4). It was discovered that all the surfaces of the materials (silica glass, aluminum alloy, 
mortar and wood) investigated far from crack initiation point, are self -affine. In addition, it was observed that the Hurst 
exponent measured along the crack direction is found to be different from the one measured along the propagation direction. 
Furthermore, ductile fracture surface had the larger fractal dimension compared with the brittle type fracture surface. 
The measurement of the porosity in aluminum cast alloys using fractal analysis was done (5). They observed that fractal 
analysis can be applied to the porosity measurement to describe the shapes of the pores in the aluminum silicon cast alloys 
using two dimensionless parameters, roughness, D and sphericity, (3. 

Spatial data analysis was used to study the distribution of micro porosity in 7050-T745 1 aluminum plate alloys (6) 
and it was observed that the micro porosity was not randomly distributed throughout the sample. Similarly, another study on 
the distribution of micro porosity in an A356 aluminum alloy plate casting was done (7). From the study, it was concluded 
that the porosity was not randomly distributed but clustered, and that the clustering tendency was relatively constant 
throughout the casting. The patterns observed in the porosity distribution maps shown in Figure 1 was used to categorize the 
pores into random, regular, clustered, and clustered with random background (8). In the present work, the intention is to use 
the fractal analysis to characterize the pores in Al - 20%wtSi, and Al - 20%wtCu cast aluminum alloys and to identify the 
points of crack initiation using MRS and SPP methods. 
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Fig. 1 : The four common types of spatial point patterns (a) random, (b) regular, (c) clustered, (d) clustered superimposed on 

random background. 

II. Materials and Methods 

lkg (80% by proportion) of commercial aluminum (99.7% pure, by weight) and 250g each (20% by proportion) of 
Si and Cu, obtained from the Federal Institute of Industrial Research Oshodi (FIIRO), were prepared for the casting. Due to 
closeness in melting temperatures, pot A which contains Al+Si, with its contents was lowered into the furnace and the 
mixture melted. 

In pot B, the alloying element Cu was first charged in the pot because of its high melting point of about 1108°C. 
When the melting point was attained, the base metal Aluminum with melting point of 660°C was added. 

The molten alloys were poured into the prepared mould after cooling for three minutes (to avoid splashing). The 
two resulting cast samples in rod forms were removed from the mould after three hours (to allow for effective cooling). 

Application of fractals 

Fractal geometry was developed (10). Its principle is universal in any measurement and has been previously used to 
numerically describe complex microstructures including graphite flakes and nodules (1, 11, 12, 13, and 14). The 
Mathematical basis for measuring chaotic objects with the power law modified is adopted in this work. The basic equation is 
as follows: 

P = P E S D1 Eql 

(l < D < 2 and 5 m < 8 < 5 U ) Eq 1 

Where P E is the measured perimeter, P is the true perimeter, S is the yardstick, 

From this expression, it can be deduced that the true perimeter is actually a function of the yardstick for 
measurement. The smaller the yardstick used, the more accurate the measurement. The fractal dimension, D, therefore 
describes the complexity of the contour of an object. It can be more practically called its roughness (5). 

When 5 < 5 m , the measurement is not sensitive to the yardstick chosen, therefore giving a smaller value of the slope, 
while when S > 8m, the size of the yardstick exceeds that of the individual feature being measured so that the measurement 
loses meaning because the object falls below the resolution limit of the yardstick used for measurement (1). Sphericity, /?, 
another dimensionless number, is used together with the fractal dimension, D; to describe the shape of the pores formed (5). 
It can be expressed as 

/3 = 4xA T /P 2 Eq2 

(0 < p < 1 and 1 < D < 2) 
Substituting equation (Eq 1) in (Eq 2) gives 

j3 = (4xA T /P?)s 2(l - D) Eq3 

(0 < p < 1 and 1 < D < l) 
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Where A T is the total pore area 

When /? = 1 and D = 1, a perfect circular shape is formed by the pore in the microstructure. 
As fi decreases, the shapes become more elongated showing a departure from perfect sphere. 
The locations of 1 < D < 2 represent less regular shapes. 



D Nodule-like Flake-tike 




1>- a 

Fig. 2: Illustration of development of irregular shapes based upon Euclidean circle or rectangle. All the shapes have the 

same area. Source: Lu and Hellawell, (1994) (1) 



Where a is the Shape Factor. 

In Figure 2 the more complex shapes for each individual island present increasing ranges of local 
curvature and the fractal dimension, D, increases. By definition, the shape factors also change, as they do for Euclidean 
shapes, but not in a readily calculated manner, because the shape factor (a) is now a function of D. 
Area of a pixel or yardstick = L x B 

Area of the total pore A T = Area of yardstick x Number of yardsticks 
To calculate the perimeter P of the pore, the Slit Island Method (SIM) 
(15) Introduced by Mandelbrot (1983) was used. It is expressed as: 

log gJ P = 0.5Dlog e A r Eq4 

log e P = log e A T % 



P = A* 



The Computer Program 

Using the equations (Eq l)-(Eq 4), an interactive Matlab program was developed to obtain the numerical values of 
the fractal dimension D and the sphericity /?. To develop the program the box counting method was used with a counter 
incorporated into the program and the small boxes or pixels occupied by the pore outlines are counted. In all four pixels 
(2x2pixels, 4x4pixels, 8x8pixels and 16xl6pixels) and four grid sizes (200x200, 100x100, 50x50 and 25x25) were selected. 
The selections were made for better resolution and to obtain accurate values. A flowchart showing the various stages and the 
subroutines in the computer program was drawn as shown in Fig. 3. 

MRS and SPP Methods. 

The first stage involves the division of the microstructures into four quadrants (lower left, upper left, lower right, 
and upper right) as shown in Figure 4. The second stage is the random selection of two pores from each quadrant while the 
third stage is the purposive selection (purposive sampling) of the "worst" and the "best" pores from the eight pores selected 
from each microstructure. The fourth stage is the categorization of the porosity distribution map into random, regular, 
clustered, and clustered with random background. Fifth stage is the discrimination between the shrinkage and the gaseous 
pores. In this stage, the patterns described in stage four are associated with different types of porosity. 
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Change colour 
of image to 
black 




Compute 
XdogtO(x) and 
Y=log10(count) 



Display the 
graph of x 
it Y 



Fit line to plot 
and calculate the 
fractal dimension 



(2) 
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(2) 



Display 
figure with 
200-X-200 
grid 



X1=0,2,4....400 
Y1=0,2,4....400 



Compute Y=X multiply by Array 1 by 201 of ones. Hold the 
current graph in the figure and plot the graph of X1 against Y. 
Also compute X=Y multiply by Array 1 by 201 of ones. Hold the 
current graph in the figure and plot the graph of X against Y1 



Display 
figure with 
100-X-100 
grid 



X1=0,4,8....400 
Y1=0,4,8....400 



Compute Y=X multiply by Array 1 by 101 of ones. Hold the 
current graph in the figure and plot the graph of X1 against Y. 
Also compute X=Y multiply by Array 1 by 101 of ones. Hold the 
current graph in the figure and plot the graph of X against Y1 



(3) 



Compute Y=X multiply by Array 1 by 51 of ones. Hold the 
current graph In the figure and plot the graph of X1 against Y. 
Also compute X=Y multiply by Array 1 by 51 of ones. Hold the 
current graph In the figure and plot the graph of X against Y1 



Display 
figure with 
25-X-25 
grid 



X1=0,16,32....400 
Y1=0,16,32....400 



Display 
figure with 
50-X-50 
grid 



Compute Y=X multiply by Array 1 by 26 of ones. Hold the 
current graph in the figure and plot the graph of X1 against Y. 
Also compute X=Y multiply by Array 1 by 26 of ones. Hold the 
current graph in the figure and plot the graph of X against Y1 



X1=0,8,16....400 
Y1=0,8,16....400 



End 



Fig. 3: Flow chart of the computer program 
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Upper left 


Upper right 


Lower left 


Lower right 



Fig. 4: The Multi-Stage Random Sampling Method of Dividing a Microstructure into four Quadrants. 

III. Results and Discussions 

Figure 5 shows the microstructures of the two alloy samples used in this work and the isolation of the pores in these 
samples. The shapes of the pores observed in the microstructures are summarized in Table 1 while the samples of the pores 
in the microstructures are shown in Figure 6. The pores are all-irregular in shape and are either nodule-like or flake-like. 

It was observed that, for the as-cast samples, the predominant pores in Al-20%wtCu are the nodule-like pores while 
the flake-like pores dominate in Al-20%wtSi (see Table 1). The formation of the nodule-like or flake-like shape is due to the 
nucleation and growth kinetics usually compounded by composition variations and the concentration of elements such as 
copper, zinc, magnesium, and silicon in the samples. 

Table 1: Sha pes of the pores and the dominating shapes in the microstructures of as- cast samples 



Cast Samples 


Shapes of the Pores 




Major Pores 


Minor Pores 


Al-20%wtCu 


Nodule-Like 


Flake-like 


Al-20%wtSi 


Flake-like 


Nodule-Like 




a: Al-20% wtCu b:Isolation of the pores in Al-20%Cu 
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c: Al-20%wtSi 



Legend: 
Dark spot . . . 
Grey spot . . . 
White spot. . . . 



Fig. 5: 



pore 

intermetallic particle 
Al-matrix 



d: Isolation of the pores in Al-20%Si 

Microstructures of as-cast aluminum alloys 




ml 




Fig. 6: Samples of pores from the microstructures 



PERCENTAGE POROSITY 

The percentage porosity, ratio of the pore area to the total area, in each of the as-cast sample, is illustrated in Table 
2. It was observed that for the as-cast samples, Al-20%wtSi has the highest number with percentage porosity of 9.92% while 
Al-20%wtCu has the lowest number with percentage porosity of 4.48%. This is because of the different composition and 
concentration of the alloying elements. The large atomic radius of Copper compared to that of Silicon ease the formation of 
pores in the interstitial spaces created in Al-20%wtSi. Copper has atomic radius of 1.57 A, while Silicon has 1.46A. 



Table 2: Percentage Porosity 
As-Cast 



Al-20%wtCu 



4.48% 



Al-20%wtSi 



9.92% 



Fractal analysis, MRS and SPP Methods 

Table 3 presents the fractal dimension, sphericity and pore location for as-cast Al-20%wtCu alloy. The porosity 
distribution map (Figure 7) gives the "best" of the pores observed as the pore with D =1.0028 and P = 0.0459 while the 
"worst" of the pores is that with D =1.0678 and P = 0.0135 corresponding to pores numbers N8 and N7 respectively. Using 
spatial point data analysis, the pores in Figure 7 have regular spatial point pattern and the porosity distribution map 
represents gas porosity because gas pores are found at a distance from their immediate neighbours due to depletion of 
hydrogen gas in the area surrounding each pore. 
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Table 3: Values of fractal dimension D, and sphericity □, for (AL20%wtCu) As- Cast Alloys 



S/n 


Alloys 


Fractal Dimension D 


Sphericity □ 


Pore Location 


Nl 


Al-20%wtCul 


1.0462 


0.0242 


Upper right 


N2 


Al-20%wtCu2 


1.0218 


0.0266 


Upper right 


N3 


Al-20%wtCu3 


1.0966 


0.0827 


Upper left 


N4 


Al-20%wtCu4 


1.0346 


0.0228 


Upper left 


N5 


Al-20%wtCu5 


1.0080 


0.0195 


Lower left 


N6 


Al-20%wtCu6 


1.0574 


0.0250 


Lower left 


N7 


Al-20%wtCu7 


1.0676 


0.0135 


Lower right 


N8 


Al-20%wtCu8 


1.0028 


0.0459 


Lower right 



It can be deduced from Table 4 (for Al-20%wtSi alloy as-cast) and the porosity distribution map (Figure 8), that the 
"best" of the pores observed is the pore with D =1.0847 and P = 0.0579 while the "worst" of the pores is that with D = 
1.2846 and P = 0.0011 corresponding to pores numbers N6 and N2 respectively. Figure 8 also has pores with random spatial 
point pattern because of the closeness of the pores to their immediate neighbours. The porosity distribution map therefore 
represents gas porosity. The implication of having only regular and random spatial point patterns, representing gas porosity 
in as-cast samples, is that the pores cannot easily link the nearest neighbour therefore making crack initiation difficult. 
Another implication is that if failure of the material will occur it will start at the locations with the worst pore shapes. 
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£ 
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0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 

Sphericity 

Fig. li Porosity distribution map for AI-20%wtCu As-cast Sample 



Table 4: Values of fractal dimension D, and sphericity □ , for (A120% wtSi) As-cast Alloys 


S/N 


Alloys 


Fractal Dimension D 


Sphericity □ 


Pore Location 


Nl 


Al-20%wtSil 


1.0510 


0.0348 


Upper right 


N2 


Al-20%wtSi2 


1.2846 


0.0011 


Upper right 


N3 


Al-20%wtSi3 


1.1087 


0.0074 


Upper left 


N4 


Al-20%wtSi4 


1.1278 


0.0104 


Upper left 


N5 


Al-20%wtSi5 


1.1633 


0.0057 


Lower left 


N6 


Al-20%wtSi6 


1.0847 


0.0579 


Lower right 


N7 


Al-20%wtSi7 


1.0440 


0.0114 


Lower right 


N8 


Al-20%wtSi8 


1.1253 


0.0126 


Lower left 
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Fig.S: Porosity Distribution Map for AI-20%TvtSi As-cast Sample 

IV. Conclusions 

♦♦♦ The noodle-like pores are the major pores in as-cast Al-20%wtCu while the flake-like pores are the major pores in Al- 
20%wtSi. 

♦♦♦ The MRS and SPP methods revealed that if crack initiation will occur it will start in: 

• Lower right region of as-cast Al-20%wtCu because it contains the "worst" pore shape with D = 1.0676 and P = 0.0135. 

• Upper right region of as-cast Al-20%wtSi because it contains the "worst" pore shape with D = 1.2846 and P = 0.001 1. 
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